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Abstract — In this work, an ensemble of economic interacting 
agents is considered. The agents are arranged in a linear array 
where only local couplings are allowed. The deterministic dynam- 
ics of each agent is given by a map. This map is expressed by two 
factors. The first one is a linear term that models the expansion of 
the agent's economy and that is controlled by the growth capacity 
parameter. The second one is an inhibition exponential term that 
is regulated by the local environmental pressure. Depending on the 
parameter setting, the system can display Pareto or Boltzmann- 
Gibbs behavior in the asymptotic dynamical regime. The regions 
of parameter space where the system exhibits one of these two 
statistical behaviors are delimited. Other properties of the system, 
such as the mean wealth, the standard deviation and the Gini 
coefficient, are also calculated. 

I. Introduction 

Nowadays it is well established that the wealth distribution 
in western societies presents essentially two phases. This 
means that the whole society can be split in two disjoint parts 
in which the richness distribution is different in each of them. 
Thus, Dragulescu and Yakovenko [1] have found, over real 
economic data from UK and USA societies, that one phase 
presents an exponential (Boltzmann-Gibbs — > BG) distribution 
which covers about 95% of individuals, those with low and 
medium incomes, and that the other phase, which is integrated 
by the high incomes, i.e., the 5% of individuals, shows a 
power (Pareto) law. Different dynamical mechanisms for the 
interaction among agents in multi-agent economic models have 
been proposed in the literature in order to reproduce these two 
types of statistical behavior. 

In the Dragulescu and Yakovenko [2] model, a set of N 
economic agents exchange their own amount of money itj, 
with i = 1, 2, • • ■ ,N, under random binary interactions, 
(ui,Uj) — > (u^u'j), by the following exchange rule: 

u\ = e(ui + Uj), (1) 
= e{ui + Uj), (2) 

with e = (1 — e), and e a random number in the interval (0, 1). 
Starting the system from an initial state with equity among 
all the agents, it evolves towards an asymptotic dynamical 
state where the richness follows a BG distribution [3]. This 
model is equivalent to the saving propensity model introduced 



by Chakraborti and Chakrabarti [4] for the particular case in 
which the agents do not save any fraction of money before 
carrying out each transaction. A similar model was introduced 
by Angle [5]. In this case, the random pair of economic 
agents, (ui,Uj), exchanges money under the rule: 

u\ = Ui - Au, (3) 

u' } = uj + Au, (4) 

where 

Au = etoui, (5) 

with e a random number in the interval (0,1). The exchange 
parameter, oj, with < oj < 1, represents the maximum 
fraction of wealth lost by the agent Ui in the transaction. When 
oj = 0.75, the exponential distribution of the wealth is also 
found for this model. 

In order to recover a power law behavior, these models can 
be reformulated in a nonhomogeneous way. These are models 
in which each agent i can have assigned by some random 
function a different value oji of the exchange parameter oj in 
the Angle model [5] or a different saving propensity in the 
Chakraborti and Chakrabarti model [4]. Under quite general 
conditions, robust Pareto distributions have been found at large 
u-values in these reformulations [4], [5]. 

Randomness is an essential ingredient in all the former 
models. Thus, agents interact by pairs chosen at random, and 
these pairs exchange a random quantity of money in each 
transaction. Moreover, the transition from the BG to the Pareto 
behavior requires the change of the structural properties of 
the system. It can be reached, for instance, by introducing a 
strong inhomogeneity that breaks the initial indistinguishabil- 
ity among the units of the ensemble. From a practical point of 
view, this could seem an unrealistic approach since we need 
to conform very different setups in order to mimic the two 
statistical behaviors. On the other hand, interactions among the 
economic or social agents of a real collectivity are not fully 
random. In fact, they are driving in the majority of transactions 
by some kind of mutual interest or rational forces. Hence, it 
would be useful to dispose of a multi-agent model with the 
ability of displaying BG and Pareto behaviors emerging from 
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an asymptotic dynamics where determinism would play some 
role in the evolution of the system. Now, we proceed to show 
with some detail the properties of one of these models [6] 
recently introduced in the literature. 

II. The Model 

The model [6] consists of a linear array of N interacting 
agents with periodic boundary conditions. Each agent, which 
can represent a company, country or other economic entity, is 
identified by an index i, with i = 1 • ■ ■ N and N being the 
system size. Its actual state is characterized by a real number, 
Xi, denoting the strength, wealth or richness of the agent, with 
Xi <E [0, oo). The system evolves in time synchronously, and 
only interactions among nearest-neighbors are allowed. Thus, 
the state of the agent, at time i+1 is given by the product 
of two terms at the precedent time i; the natural growth of the 
agent, Tix\, with its own growth capacity parameter, ri, and 
a control term that limits this growth with respect to the local 
field 'J* = \{x\_ x + x\ +l ) through a negative exponential 
function with local environmental pressure af. 

xj +1 =nx\ exp(- \x\- ai ^\ \). (6) 

The parameter r,; represents the capacity of the agent to 
become richer and the parameter describes some kind of 
local pressure [7] that saturates its exponential growth. This 
means that the largest possibility of growth for the agent is 
obtained when Xi ~ a^', i.e., when the agent has reached 
some kind of adaptation to the local environment. In this note, 
for the sake of simplicity, we concentrate our interest in a 
homogeneous system with a constant capacity r and a constant 
selection pressure a for the whole array of sites. 

If all the agents start with the same wealth, the index i can 
be omitted, x\ = x l and "f* = x l , and the global evolution 
reduces to the following map, 

x t+1 = nc*exp(- | 1 - a \ x*). (7) 

The above map can be easily analyzed by standard tech- 
niques. In fact, the parameter a could be removed by doing the 
change of variable y* =| 1 — a \ Xt, and obtaining the generic 
map y t+1 = ry l exp(— y*). For r < 1 the system relaxes to 
zero and for r > 1 the dynamics can be self-sustained deriving 
toward different types of attractors. It displays all kind of 
bifurcations known for this type of maps [8], except, evidently, 
for the singular case a = 1. For instance, when r > 1 the fixed 
point is xq = logr/ | 1 — a |. This point becomes unstable by 
a flip bifurcation for r = e 2 . For increasing r the whole period 
doubling cascade and other complex dynamical behavior are 
obtained. However, it can be shown that such evolving uniform 
states are unstable. When a perturbation is introduced in the 
initial uniform state or, in general, when the initial condition 
is a completely random one, the asymptotic dynamical state 
of the system is found to be more complex. 



III. Some Results 

We perform the statistical study of the system in the 
parameter space (a, r). For all the simulations, the system 
size is N = 10 5 and the initial conditions in the array are 
random values in the interval (1, 100). Also, a transient of 
10 4 iterations is completed before arriving to the asymptotic 
state where all the measurements are done. At this point, if 
necessary, the average is done over the next 100 iterations 
after the transient, and this result is newly averaged with the 
same process over 100 different realizations of the system. 
Following this method, different statistical quantities, such 
as the mean wealth, the standard deviation and the Gini 
coefficient, have been calculated on the model, and they will 
be shown in the communication to be presented in NOMA'07. 

As an example, we advance here the possibility of existence 
of the two statistical behaviors in which we are interested in. 
Hence, if the number of agents, P(x), with a wealth x is 
represented in a semi-log plot, a BG distribution is found for 
the values a = 0.6 and r = 4 (Fig. la, where values of maps 
have been taken at time t — 10 4 ). The exponent /i of this 
distribution, P(x) ~ e~^ x , is fi = 0.26. A thermodynamical 
simile can be done for this type of behavior by defining a 
kind of 'temperature', h = 1/fi, that is related with the mean 
wealth of each agent in the ensemble. In this particular case, 
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Fig. 1. (la): Semi-log plot of the number of agents, P(x), with a wealth x 
for a = 0.6 and r = 4 in the asymptotic regime at time t = 10 4 . The scaling 
of this exponential behavior is 0.26. (lb): Log-log plot of P(x) for a = 0.92 
and r = 8 at time t = 10 4 . The exponent of this power law behavior is 2.84. 

Also, a power law behavior in the statistical behavior of 
system (|6]l is possible, for instance, for a — 0.92 and r = 
8. If the number of agents, P(x), with a wealth x is now 
represented in a log-log plot, a Pareto distribution is unmasked 
[6] (Fig. lb, where values of maps have also been taken at time 
t = 10 4 ). The exponent a of this distribution, P(x) ~ x~ a , 
is in this case a = 2.84, which is in clear agreement with the 
exponents derived from real economic data [1], [9], [10]. Thus, 
the exponent a found in the cumulative probability distribution 
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of incomes, whose scaling behaves as x~ a ~ x~ a+ , is about 
2.1 for the UK and 1.7 for the US economies [1]. Pare to 
himself proposed a value of 1.5, Levy and Solomon [9] found 
a value of 1.36 for the distribution of wealth in the Forbes 400 
and Souma [10] found 2.05 for the high income distribution 
in Japan. 

IV. Conclusion 

Let us conclude this note by remarking that our present 
model © shows, at least, two interesting statistical behaviors, 
those of Boltzmann-Gibbs and Pareto types, in its asymptotic 
state. The procedure to get these two regimes does not require 
of any structural change in the system. Only, it is necessary to 
tune some adequate pair of values of the external parameters 
(a, r) for which the system displays one of those types of 
behavior. This property could be an advantage respect to other 
models in the literature that need to perform a strong change in 
the conformation of the initial ensemble in order to present one 
of the two behaviors, exponential or power law. Other relevant 
property of model (O is its complete determinism. There is 
no any kind of random ingredient in its evolution. Of course, 
this does not forbid in any way that the asymptotic state of the 
system can show some degree of spatio-temporal complexity, 
including very large fluctuations of the agent's wealth. As 
indicated by Yakovenko in his last review on this subject 
[11], let us finish by saying that it seems well proved that 
western societies consist of two different classes characterized 
by different distribution functions. However, the most part of 
theoretical models on this subject do not produce two classes, 
although they do produce broad distributions. We hope that 
the model here studied can be used in the next future as a tool 
in the agent-based theory capable of simulating two economic 
classes populations. 
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